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Abstract 

The statistical production of antibaryons is considered within the canonical ensemble 
formulation. We demonstrate that the antibaryon suppression in small systems 
due to the exact baryon number conservation is rather different in the baryon-free 
(B = 0) and baryon-rich (B > 2) systems. At constant values of temperature and 
baryon density in the baryon-rich systems the density of the produced antibaryons 
is only weakly dependent on the size of the system. For realistic hadronization 
conditions this dependence appears to be close to B/{B + 1) which is in agreement 
with the preliminary data of the NA49 Collaboration for the p/n ratio in nucleus- 
nucleus collisions at the CERN SPS energies. However, a consistent picture of 
antibaryon production within the statistical hadronization model has not yet been 
achieved. This is because the condition of constant hadronization temperature in 
the baryon-free systems leads to a contradiction with the data on the p/ir ratio in 
e"'"+e~ interactions. 
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1. Introduction 

Among the first models of multiparticle production in liigli energy interactions were 
statistical models 0, In the last decade a significant development of these models 
and the extension of the area of their applicability took place. The main reason for 
this is a surprising success of the statistical approach in reproducing new experimental 
data on hadron multiplicities in nuclear (A+A) [^] and elementary (e"'"+e~, p+p, p+p) 
collisions P). One of the important results of the analysis of hadron yield systematics 
at high energies (SPS and higher) done within the statistical models is the approximate 
independence of the temperature parameter T = 160 -v- 190 MeV from the system size 
and collision energy [0]. This result can be attributed to the statistical character of the 
hadronization process. 

The statistical models are based on the key assumption that all microscopic states of 
the system allowed by the conservation laws are equally probable. Calculations within 
these models are straightforward when the mean number of particles of interest is large and 
consequently it is enough to fulfill the conservation laws in the average sense, i.e., for the 
macroscopic state. This is achieved by the introduction of parameters: the temperature T 
and chemical potentials /Xj, which control correspondingly the average values of the system 
energy density and its material content (e.g., baryon number, strangeness and electric 
charge). In this case (grand canonical ensemble, g.c.e.) the mean particle multiplicities 
are just proportional to the volume V of the system. The particle density and the ratio 
of the multiplicities of two different particles often used for the comparison with the data 
are volume independent. 

This simple volume dependence is however not valid any more for a small system in 
which the mean particle multiplicity is low. In this case (canonical ensemble, c.e.) the 
material conservation laws should be imposed on each microscopic state of the system. 
This condition introduces a significant correlation between particles which carry conserved 
charges. The correlation reduces the effective number of degrees of freedom and conse- 
quently leads to the c.e. suppression of the 'charged' particle multiplicity when compared 
with the result of the calculations done within g.c.e. A 'neutral' particle (a particle which 
does not carry conserved charges) does not feel the c.e. suppression in the small system, its 
mean multiplicity remains proportional to the volume. Therefore, a different dependence 
on the system volume is expected within statistical models for 'charged' and 'neutral' par- 
ticles. The magnitude of the c.e. suppression of the 'charged'- 'anticharged' pair creation 
increases with the mass of the lightest hadron needed to compensate the particle charges. 
Therefore, a strong c.e. suppression may be expected for antibaryon production as the 
mass of the lightest baryon is m = 938 MeV which is much larger than the value of the 
temperature parameter found in the hadronization models. The c.e. suppression is still 
rather essential for strange particle production (see, e.g., p|]). 

The above expectation seems to be violated by the preliminary data on antiproton 
production presented recently by the NA49 Collaboration [Q. The p/vr ratio is found 
to be approximately the same for p+p interactions and central Pb+Pb collisions at 158 
A-GeV [antiproton scaling). Thus the hadronization volume increases but the effect of 
the c.e. suppression is not observed. The c.e. suppression can be expected in p+p inter- 
actions because of the key difference between antiproton and pion: the antiproton carries 
baryon number in addition to the electric charge carried by both particles. In order to 
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compensate the electric charge of a produced particle it is enough to create an additional 
charged pion. As the pion mass is much smaller than the nucleon mass, a significantly 
stronger c.e. suppression is expected for antiproton production than for the pion pro- 
duction, consequently it should lead to a strong violation of the experimentally observed 
scaling. Thus the crucial question is whether the antiproton scaling can be understood 
within the statistical model of hadron production in which the condition of exact baryon 
number conservation is imposed. We note that the antiproton multiplicities in high energy 
collisions were shown to approximately agree with the predictions of statistical models 
10, H, P] . However different versions of the models with different parameters were used to 
fit various sets of data. Thus the question whether consistent description of the antiproton 
data in the statistical model is possible is still opened. 

The importance of the exact treatment of the material conservation laws within sta- 
tistical models of strong interactions was first pointed out by Hagedorn [|l^ (see also 
Refs. []TT], [T^). Subsequently a complete treatment has been developed (see, e.g., and 
references therein) and applied to analyze the hadron yields in elementary collisions ^ . 
In this letter we derive explicit analytical formulae to study the role of the exact material 
conservation laws within the statistical model of hadronization and we use them to discuss 
the antiproton scaling ohseTved experimentally in Pb+Pb collisions at 158 A-GeV 0. We 
also discuss the data on the p/ir ratio in e"'"-|-e~ interactions |]14|, |1^ within the statistical 
hadronization model. 



2. Model formulation 

Let us consider the system of baryons '6' and antibaryons 'a' with total baryon number 
B as the Boltzmann ideal gas in the volume V, at temperature T. The c.e. partition 
function is 



Z{T,V,B) 



n 



E - E - 



(i)i 



where the index j runs over all (non-strange) baryon states N, A, N* , and the single 
baryon (antibaryon) partition function reads 



Zj = Zj{T,V) 



(2vr) 
27r2 



(2) 



T m] K,{m,/T) ^ V f,{T) . 



The baryon mass and the baryon degeneracy factor are denoted here by mj and gj, 
respectively. Auxiliary parameters A^"''' and A^-'-' are introduced in order to calculate the 
mean number of baryons and antibaryons and they are set to unity in the final formulae. 
By expressing 6k as 

1 /-Ztt 

5K(n) 



2tt Jo 



-incj) 
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Eq. (|T]) becomes 
Z{T,V,B) 



1 r^^ 
1 r^^ 

271 Jo 



-iB(j> 



n E E 



(A^^'^^, e'^)^^" {Xi^hj e-^<^)^' 



e-*^^ exp 



(3) 



This form of the c.e. partition function allows one to derive the mean numbers of baryons 
and antibaryons 



(i) 



b ' Ab=Aa=l 

'd\ogZ\ 



(i) 



Zi 



Z{T,V,B-1) 
Z{T,V,B) 

Z{T, V,B + 1] 
Z{T,V,B) 



(4) 



(5) 



For Xb = Xa = I the partition function can be presented as the modified Bessel function 

1 r'^T^ 



Z{T,V,B) = — r d(t)e-'^^ exp(2z cos( 
27r JO 



Ib{2z) 



(6) 



where z = This yields final expressions for the mean number of baryons and 

antibaryons 



^ Ib-i{2z) 
Ib{2z) 



{2z) 



Ib{2z) 



(7) 



As the exact baryon number conservation is imposed on each microscopic state it is 
evidently fulfilled also by the average values (^: 

{N,) - w = 5:(ivi^-)) - 5:(iv(^-)) = 5 , (8) 

j 3 

as indeed can be easily seen from the identity — = 2nln{x)/x [|I^. Eq. (|^) 

is valid for all combinations of B and z values. For a specific case of 5 = in the nucleon- 
antinucleon gas (i.e., no resonances included) our results (|^ are reduced to the result of 
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The c.e. expressions for the mean number of baryons and antibaryons can be further 
simplified for the two hmiting cases: z <^ 1 (small systems) and z 1 (large systems). 
Using the representation of /„ as the infinite series |T^ 



Wz) 



z 



n+2k 



one obtains for small systems 



2j £j_ _l_ ' ^ 



- + oizrz'), (iVP) = + oiz,.z'). 



B + 



(9) 



4 



The dependence (No) oc V"^/ {B + 1) is therefore observed from Eq. (|^) for the antibaryon 
yield in small systems. Such a dependence can be intuitively understood from the ki- 
netic picture of the baryon-antibaryon pair creation and annihilation. Let's consider the 
time dependence Nait) (the time averaging of Na{t) should reproduce our statistical av- 
erage (Na)). At each time moment the value of Na{t) equals or 1 (configurations with 
Nait) > 2 can be safely neglected as {Na) ^ 1). The kinetics of the evolution of Na{t) is 
defined by the frequency uj of the baryon-antibaryon pair creation and the life-time At 
of the produced baryon-antibaryon pair: {Na{t)) = uAt. As a baryon-antibaryon pair is 
locally produced in any point of the system we have u (x V. Being produced as baryon- 
antibaryon pair the antibaryon can be then locally annihilated by any baryon existing in 
the system. Therefore, At oc I/ub {ub is the density of baryons). hb = ^-/V loi B = Q 
(as only one baryon exists in the volume V^), and Hb = (B + l) /V if B baryons are present 
in the system before the baryon-antibaryon pair creation. These lead to the dependence 
of the antibaryon number on V and B as given by Eq. (^. 

For large systems {z ^ 1) the c.e. becomes equivalent to the g.c.e. where the partition 
function and the average number of baryons and antibaryons are calculated as 

Z{V,T,fiB) = E ^^P(^) Z{T,V,b) = exp(ze^-/^ + ze"'^-/^) , (10) 

(Ar(^)) = zj exp(/i5/T) , (Ar(^)) = zj exp(-/.B/T) . (11) 
Here /ib is a baryon chemical potential which using Eq. (|^) is defined for B = {b) >0 as: 




expifiB/T) = _ + Wl+ _ . (12) 



Note that the function fj = fj{T) introduced in Eq. (|^) has the physical meaning of the 
density of the j-th baryon and antibaryon in the g.c.e. formulation for fiB = 0. Using the 



uniform asymptotic expansion of the modified Bessel functions at n ^ cxo |[T^ 
1 exp{nri) 



(l + a;2)i/4 



l + o( - 



n 



. r^^Vl+^ + ln , (13) 

1 -I- V 1 + a;^ 

results (pT|) and (|1^) for the g.c.e. are also easy to obtain from the c.e. (|^) in the ther- 
modynamical limit V oo, B ^ oc with B/V = pb = const{V). 

In order to remove a 'trivial' linear dependence of the particle multiplicities on the 
system volume it is convenient to make a comparison between the particle ratios from the 
model and the experimental data. In the statistical model the multiplicity of any 'neutral' 
meson state M is just proportional to the volume {{Nm) = zm = V fuiT)) for both small 
(c.e.) and large {g.c.e.) systems. Therefore, the system volume dependence of the ratio 
{Na) I {Nu) at fixed temperature is the same as for the antibaryon density. From Eqs. (|^) 
and (|Tl|) the antibaryon densities for the c.e. and g.c.e. are equal to 

k^^^) , _ f,) , . ^B+i(2z) , Ib+i{xB) ^ ^ X B 



n 



V '^^ " '^^ Ib{2z) Ib{xB) ~^'2B+1' ^ ' 



n 



(J) 



X 



fj exp{-pB/T) = fj ^ , , , (15) 



y \gce — '"a \gce jj ---i-v t^o/ - / ]^ _j_ y/Y+ 



X 
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where / = J2j fj, ^ = '^z/B = 2f /pb, and the last approximation in Eq. (p^ ) is vahd for 
small systems only. Note that introducing the variable x we have transformed the finite 
size l^-dependence of the ce. density ([T^) into its dependence on the baryon number 
B. Eqs. give us the primary thermal density for all individual antibaryon states 

j. Each non-strange resonance (anti)baryon state decays finally into (anti)nucleon plus 
meson(s). Therefore, the total (primary plus resonance decay) antinucleon density equals 
to the total thermal antibaryon density, Ua = Y,j n-a^ and is given by Eqs. ([l4| , p!5D with 
the substitution of fj by a sum / = J2j fj- 

For the purpose of the following discussion we define a canonical suppression factor 

It quantifies the antinucleon suppression due to the exact baryon number conservation. 
We note also that the suppression factor Fes 0T6| ) is the same for any individual antibaryon 
state. 

3. Discussion 

The results derived in the previous section are used here to discuss antibaryon pro- 
duction in high energy collisions. 

In the B = case the baryon and antibaryon densities are equal and Eqs. (|14D and 
(H) yield 

«-U^4"U^/,^-/,. «-U^ 4"U^A, (17) 

where the approximation for the ce. density is valid for small system only. The canonical 
suppression factor (|l^) for B = is equal to 

FO = = fV . (18) 

The approximation in Eq. (^) is valid for small system only {z = fV <^ 1). 

The behavior of the canonical suppression factor F^^ (|T8|) is shown by the solid lines 
in Fig. 1 for T = 160 MeV, 170 MeV and 180 MeV, assuming that / is the sum of fj over 
all non-strange baryons. The lines start from V = 5 fm^, which is approximately equal 
to the estimate of the hadronization volume for e"*" -|- e~ interactions at ^/s = 29 GeV 
0. One observes (see Fig. 1) a strong ce. suppression of the (anti)baryon density. For 
T = 160 MeV the (anti)baryon density increases by a factor of 10 from its value at 

= 5 fm^ to its V — >■ oo g.c.e. limit. For the small systems the (anti)baryon density 
increases approximately linearly with V, i.e., the (anti)baryon multiplicity for the small 
systems is proportional to V"^. The ce. suppression becomes less pronounced and the 
volume region with linear increase of the (anti)baryon density is reduced for increasing 
temperature. 

Let us now turn to the antibaryon production in baryon rich system. In the analysis of 
data on particle multiplicities in p-|-p, p-|-A and A-l-A collisions one usually assumes that 
all participating nucleons in the collisions (wounded nucleons) take part in the statistical 



6 



hadronization of the system. It means that in the analysis of the NA49 results on antipro- 
tons from p+p interactions to central Pb+Pb collisions at 158 A-GeV we should study 
statistical systems with 2 < B < 400. It was found that the mean multiplicity of pions 
per wounded nucleon increases (at the SPS collision energies) only by about 20% when 
going from p+p interactions to central Pb+Pb collisions [|1^. The pion to baryon ratio 
in the statistical model is determined by two parameters: the temperature and baryon 
density. Thus as the temperature is found to be constant (T = 175 ± 15 MeV) we con- 
clude that the baryon density at hadronization in nuclear collisions at 158 A-GeV is also 
approximately constant. 

Therefore, for the comparison with the NA49 results we study the evolution of the 
antibaryon density with increasing net baryon number B aX T = const and ps = const. 
The c.e. suppression factor (0) is found at these conditions from Eqs. (^[15]) 



_ l + y/TT^ Ib+i{xB) 2/ 

Ib{xB) ' ^ 

Its i?-dependence is plotted in Fig. 2 for several different values of the parameter x. Note 
that our assumption T = const and ps = const for statistical hadronization at different 
values of B can be substituted by a weaker one, x = const. From Fig. 2 one observes that 
the c.e. suppression of antibaryon density becomes stronger at high baryon density (i.e., 
small x). For x < 1 the c.e. suppression (|T9|) becomes close to its x — limit: 

Ffs = ■ (20) 

B + 1 



Eq. (pOD shows that the strongest c.e. suppression of the antibaryon density is for the 
B = 2 (nucleon-nucleon interactions) case and it leads to the suppression factor of 2/3. 
This moderate effect of c.e. suppression is in strong contrast with the large c.e. suppression 
(i.e., <^ 1) in the baryon-free system. A mathematical reason of this very different 
behavior for B = and B > 2 (with ps = const {V)) is due to the fact that in the latter 
case both the order of the modified Bessel functions and their arguments are dependent 
on B (i.e., on V) whereas in the B = case only the argument increases with V. 

The presence of non-zero baryon number B > has a twofold effect on antibaryon 
production. First, it suppresses the production of antibaryons: the additional factors 
exp(-/iB/T) = x/(l + VI + a;2) < 1 and 1/(5 + 1) < 1 appear respectively in the 'large' 
and 'small' systems for the antibaryon density in comparison with the B = case. On 
the other hand, the c.e. suppression effect due to the exact baryon number conservation 
becomes smaller: at fixed T and V the following inequality is always valid, F^ > Fj^^. For 
fixed B > the c.e. suppression of antibaryons becomes smaller when ps decreases and 
it disappears completely (i.e., F^ 1) in the limit Pb (and respectively V — > oo in 
order to keep the B value fixed). This is because the total number of baryon-antibaryon 
pairs becomes large due to large V. Note that in this case the last approximation in 
Eq. (|T^) is no more valid. Instead one should use the large argument asymptotic of the 
modified Bessel functions. 

Thus for i? > 2 systems at constant x = 2f / ps the c.e. suppression factor ([TP|) 
ranges between 2/3 and 1 for a; ^ 1 and between (1 — l/4a;) and 1 for x^ 1. 
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Previous analyses of hadron production at the CERN SPS indicate large baryon den- 
sities at hadronization. The typical values of T = 170 MeV and fis — 250 MeV found for 
Pb+Pb collisions |Q lead to the estimate x = 2f / pB = sh~^{fiB/T) ~ 0.5. As seen from 
Fig. 2 the c.e. suppression factor (p^Of ) for this 'small' value of x is close to its limiting 
pattern B/{B + 1) In Fig. 3 the NA49 results on the p/vr ratio in p+p and Pb+Pb 
collisions at 158 A-GeV are compared with this limiting pattern. From this comparison 
we conclude that the model of statistical production of antiprotons at hadronization in 
baryon-rich system correctly reproduces the observed antiproton scaling. 

Let us return again to the case of the baryon-free system. The statistical model cal- 
culations for e++e~ [|| and p+p interactions include large c.e. suppression effects. As 
discussed in the introduction we assume that the hadronization temperature reflects a 
universal property of the hadronization process and therefore should be collision energy 
independent. In the case of e"'"+e~ ^ interactions the hadronization volume is small and 
therefore one expects approximate proportionality to V"^ of the multiplicity of nucleon- 
antinucleon pairs but only a linear increase with V of the pion multiplicity. Therefore, 
the p/tt ratio ratio calculated within the model increases linearly with increasing pion 
multiplicity. However experimental data |jl^, contradict this expectation of the statis- 
tical model. The p/ir ratio which is plotted in Fig. 4 as a function of pion multiplicity for 
e"'"+e~ interactions at different energies, a/s = 14 + 91 GeV, is approximately constant. 

Within the discussed statistical hadronization model one can try to solve the problem 
by assuming an increase of the temperature T with decreasing volume V. The function 
/(T) strongly increases with T which allows to compensate the c.e. suppression effect (to 
keep the (anti)nucleon density, Ua = p(T)V, constant) for moderate (of about 10 MeV) 
changes of T. This indeed is observed in the fit results of the statistical model for the 
e"'"+e~ and p+p data: the increase of the volume is always accompanied with the decrease 
of the temperature parameter ^ . Thus one may argue that the hadronization condition 
T = const has to be substituted by a different criterion which should explain the decrease 
of the temperature with increasing size of the system in e^+e~ and p+p interactions. The 
constant energy per particle was recently discussed as a chemical freeze-out condition |l^ . 
The detailed study of this question is, however, outside of the scope of the present paper. 



We note only that the statistical production of heavy particles (e.g., J/ip mesons ||T9| ) 
is very sensitive to the temperature parameter. Their yields, therefore, can be used to 
clarify the problem. 



4. Summary 

The role of baryon number conservation in the calculations of antibaryon multiplicity 
within statistical model of hadronization was investigated. We derived explicit analyt- 
ical formulae for the antibaryon multiplicity in baryon-free and baryon-rich small and 
large systems. This formalism was further used to discuss antiproton scaling observed 
experimentally in A+A collisions. The statistical model with constant hadronization 
temperature correctly reproduces the weak dependence of the p/vr ratio on the system 
size in p+p and nuclear collisions at the CERN SPS energy. A description of the ratio of 
J/ip mesons to pions within the statistical hadronization model requires also a constant 
temperature parameter in p+p and A+A collisions at the CERN SPS. However, the same 
model with T = const does not give a natural explanation of the approximate indepen- 
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dence of the p/tt ratio of collision energy in e++e~ interactions. Therefore, a consistent 
description of hadron production within the statistical hadronization model has not yet 
been achieved. 
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Figure 2: The finite size -B-dependence of tlie antibaryon production in baryon ricli 
{B > 2) systems at different values of tlie variable x [x = 2z/B = 2f /pb)- Tlie solid 
lines show the c.e. suppression factor (|19D for x=l and x=b (from below to above). 
The lower dotted line corresponds to the limiting B/{B + 1) behavior (|20|). 
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Figure 3: The NA49 data on the antiproton to pion ratio {{N^) = {{N^+) + {N^-))/2) 
in p+p (square) and centrahty selected Pb+Pb (dots) coUisions at 158 A-GeV are plotted 
as a function of the mean number of wounded nucleons, (Np). The dependence of the 
ratio on (Np) expected within the statistical model, {Np) /{{Np) + 1), is shown by dotted 
line, the function is normalized to the experimental data. 
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Figure 4: The antiproton to tt" |]T5[ ratio in e~^e~ annihilation at y/s = 14, 22, 29 
35 and 91 GeV is plotted as a function of mean multiplicity of tt" mesons, {N^^o). In 
order to increase the significance of the data the antiproton multiplicity was calculated as 
average of proton and antiproton multiplicities. The dotted line indicates the proportional 
dependence of the ratio on {N^^o) expected within statistical model for small systems. The 
horizontal dashed line is shown for the reference. 
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